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Abstract 

We investigate the properties of strongly interacting heteronuclear boson-boson mixtures loaded 
in realistic optical lattices, with particular emphasis on the physics of interfaces. In particular, we 
numerically reproduce the recent experimental observation that the addition of a small fraction 
of ^^K induces a significant loss of coherence in ^''Rb, providing a simple explanation. We then 
investigate the robustness against the inhomogeneity typical of realistic experimental realizations 
of the glassy quantum emulsions recently predicted to occur in strongly interacting boson-boson 
mixtures on ideal homogeneous lattices. 
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Ultracold degenerate gases in optical lattices provide an unprecedented toolbox for realiz- 
ing experimentally what were once just toy models sketching the key features of complicated 



)bard (BH) model, orig- 
ll and later adopted for 



condensed matter systems. One prominent example is the Bose-Hu 
inally introduced as a variant of the better known Hubbard model 
the description of superfluid ^He trapped in porous media |2j]. Several years after the intro- 
duction of this simple yet challenging toy model, Jaksch and co-workers suggested that it 
could be realized in terms of ultracold bosonic gases trapped in optical lattices (sf , and were 
soon proved right by a breakthrough experiment where the hallmark superfluid-insulator 
quantum phase transition of the BH model was observed J]. 

Recently, several experimental groups directed their efforts to the realization of more 
complex generalizations of the Hubbard model, involving mixtures of particles obeying either 
the same or different statistics. Beyond their theoretical appeal, these systems are relevant 
to interesting applications such as implementation of disorder ^, ^, association of dipolar 
molecules j?], schemes for quantum computation 8] and mapping of spin arrays j^. 

So far, most of the experimental efforts on optical lattice systems have been directed to 



boson-fermion mixtures 



a 
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while fermion-fermion and boson-boson (BB) mixtures 



have been somewhat ignored. Very recently the Florence group performed an experiment 
on a harmonically trapped BB mixture of atomic ^^K and ^^Rb with strong interspecies 
repulsion [l^. Expectedly, the presence of a relevant K fraction modifies the quantum 
phase transition occurring in Rb. More surprisingly, this effect turns out to be sizeable even 
for a small overlap between the two atomic species . Strongly interacting BB mixtures are 
also the subject of a recent theoretical investigation, whose main observation is that strong 
interspecies repulsion can substitute for disorder, driving a mixture loaded in a homogeneous 



1-D lattice into metastable quantum emulsion states exhibiting glassy features 



13|- 



In the present work we introduce a unified framework for the description of lattice BB 
mixtures with strong interspecies interactions in realistic conditions and different physical 
regimes encompassing and generalizing the above- described findings [l^, 3]. In particular, 
we explain the apparently surprising observation that the coherence properties of a bosonic 
system can be reduced significantly even in the presence of a single interface with a second 
bosonic species Furthermore, we establish the range of parameters for which the intu- 
itively expected opposite behavior of increased coherence is recovered. Concerning quantum 
emulsions, we show that they are in principle compatible with the inhomogeneity arising 
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from confining potentials typical of experimental realizations, albeit in a restricted range of 
Hamiltonian parameters. Specifically, while in the homogeneous case a sufficiently strong 
interspecies repulsion ensures the occurrence of quantum emulsion states |l3|, in the exper- 
imentally relevant inhomogeneous case the difference of intraspecies repulsions turns out to 
be a fundamental critical parameter. 

The systems under concern provide a realization of the two-flavor BH Hamiltonian 



U 



V. 



(1) 



where the lattice boson operators aj- ,^-, afj, and Ufj = a^^ -afj, create, destroy and count 
atoms of type / at site j. The parameters f//, U12 quantify the intra- and inter-species 
BB (repulsive) interaction, Jf is the hopping amplitude and f/j = mjQf^j — j°)^/2 is the 
standard harmonic trapping potential felt by bosons of species / at lattice site j. By ruf, 
Qf, and jj, we denote, respectively, the mass, the trapping frequency, and the minimum 
point of the harmonic potential vjj of species /. 

Since our aim is the study of strongly interacting mixtures, it is convenient and effective 
to adopt a mean-field approach based on the assumption that the ground state of the sys- 
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tem is the product of on-site factors I^P) = Ylj li'j), \ipj) = Sm.na'^^i'^a (^'^ij j [^Ijj 1^); 
where \Q) is the vacuum state, af,j\^) = 0, and the coefficients Cn^^n2 are determined via 
energy minimization at fixed atomic populations A^^i, N2. Owing to a much lower com- 
putational demand, this mean-field approach provides qualitative results on systems that 
would be beyond the present capabilities of more quantitative numerical methods, such 
as quantum Monte Carlo (QMC), density matrix renormalization group or time-evolving 
block-decimation algorithms. 

. iamiltonians similar to that under examination have been considered previously 




same bosonic species [9 
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23|, possibly referring to two different internal states of the 



Iq to spin-1 [19| or dipolar bosons [7|, l22|. Most of the previous 



work focuses on the phase diagram of hornogeneous lattices, often adopting a mean-field 



approximation similar to ours 
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However, our approach is characterized 



by some features that have not been considered in the literature, at least simultaneously. 
First of all, our mean- field is fully site-dependent, and does not reduce to an effective single- 



site theory. This allows us to describe phase-separated systems and to consider realistic 
harmonic trapping potentials. Furthermore, we fix the bosonic populations, A^^i and A^2, 
rather than the corresponding chemical potentials, /ii and fi2- Again, this allows us to make 
direct contact with experimentally relevant situations, and avoids the "species depletion" 
problem Q, i.e. the vanishing of one atomic species from the minimum-energy state at 
fixed chemical potential. 

The method is first applied to a situation reproducing the experimental conditions in 
Ref. IJ], where a bosonic mixture of Rb and K was loaded in an optical lattice. The 
potentials trapping the two atomic species had the same kf = rufVLf but, since mi 7^ 
m2, their minima were displaced in the vertical direction: = 2g/{VL^^X)^ where g is the 
gravitational constant and A/2 is the optical lattice spacing (henceforth the subscripts 1 and 
2 will denote Rb and K, respectively). An important consequence of the interplay between 
the ensuing asymmetry and the strong interspecies repulsion is the tendency towards full 
phase separation, minimizing the number of interfaces between the two species. In fact, in 
the Florence experiment the interspecies overlap is estimated to be limited to one lattice 
site in the vertical direction. Despite the occurrence of a single phase interface, the effect 
of K on the coherence properties of Rb turns out to be sizeable [12;]. More in detail, it has 
been observed that a modest quantity of K (around 10% of Rb) reduces the coherence of Rb 
significantly, moving the superfiuid-insulator transition point to smaller values of the optical 

a naive argument based on Ref. B 
results in a prediction opposite to the observed behavior: the presence of K increases the 
local density of Rb, which would cause an increase in the coherence of the latter. 

This argument is indeed valid for most of the phase diagram of the BH model describing 
an atomic cloud loaded in a homogeneous optical lattice. However, clear exceptions are 
found in the proximity of the Mott lobes, where an increase of the (local) density -or chem- 
ical potential- results in a sharp drop in the condensate fraction. Furthermore it should 
be emphasized that such phase diagram describes a homogeneous system in the thermo- 
dynamic limit, whereas here we are dealing with a finite and inhomogeneous system. The 
site-dependent potential acts like a local chemical potential for a system with fixed total 
population. As a result, at sufficiently high ratios of interaction to kinetic energy, config- 
urations of the system can be found where superfluid and Mott-insulating domains coexist 
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251]. The density of the system assumes the so-called wedding-cake or ziggurat profile, 
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FIG. 1: (color online) a) \ai^h? vs iVi {s = 11, Ji « 0.022 f/i, ki « 7.9 • lO^^C/i); b) same as a), 
but in the presence of K (A^2 = 30, U2 ~ 0.65 Ui, J2 ^ 0.21 Ui, Uu ~ 2.22 C/i, ^2 = h c) /c of Rb 
corresponding to a) (black) and b) (gray); d) and e) visibility of Rb vs. lattice strength, f) and g) 
Configurations of 340 Rb (blue) and 30 K (red) atoms for s = 11, along with the relevant trapping 
potentials (arbitrary units). The height of each bar represents the local population Uh, whereas 
the darkness of the shading is proportional to l^j./jp. In f) we set C/12 = 0. 



the plateaus corresponding to (quasi) Mott-insulator domains. When the configuration is 
such that the topmost plateau involves a fair number of sites, the density profile responds to 
an increase in the total population according to a predictable sequence. At first a dome-like 
essentially superfiuid structure appears on top of the highest plateau. Subsequently, the 
width and height of this structure increase, leading to an increase in the system coherence. 
When the tip of the dome gets too close to the next level of the ziggurat, the dome fiattens. 



its central part turning gradual 
overall coherence of the system 



y into a plateau. Correspondingly, there is a drop in the 
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The above described single-species scenario is captured quite satisfactorily by the 



Gutzwiller mean-field approximation 



21, 



281]. We will now show that it bears a 



strict relation with the experimental observations reported in Ref. 12] about the Rb-K 
BB mixture. Fig. [T] shows results obtained from a double-species Gutzwiller mean-field 
approach where we have adopted physical parameters — J/, Uf, Ui2,^f, atomic density 



at the trap center, population ratios — in the experimentally determined range 
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29]. 
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FIG. 2: number of interfaces vs energy for A^i = = 180 particles on a 300-site ID lattice. In all 
cases Ji = J2 = O.lUi and U12 = l.bUi. 



For the sake of simplicity we have focused on a 1-D lattice as mean-field results are essen- 



tially independent of the dimensionality 30|]. Panel a) shows the local superfluid parameter 
= I P of R-b alone as a function of the relevant population A'^i and lattice 
site label, h (the darker the hue, the larger lai^/^p). The drop in the superfluid parameter at 
the trap center signals the formation of new ziggurat levels from the flattening of coherent 
domes. Panel b) shows the same quantity as in a) yet in the presence of 30 atoms of K 
(|a2,jp is not shown). The main effect of the addition of K is that the new structures of the 
(now asymmetric) ziggurat appear at smaller populations A^^i. Panel c) shows an estimate 
of the coherence of Rb measured in terms of the relevant condensate fraction /q [31] for 
the data in panels a) (black) and b) (gray). The presence of K is indeed equivalent to an 
increase in Rb population, but, given the oscillatory behavior of /q, this does not necessarily 
result in an increase of the overall coherence of Wo. A small fraction of N2 can cause either 
an increase or a decrease of /q, depending on the value of A^i. The experimental measure of 
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coherence, i.e. the so-called visibility V 



32l |. exhibits similar oscillations as in c), albeit 
with a different envelope. Panels d) and e) show the changes in V^^) produced by N2 = 30 K 
atoms, for two values of Ni. Note that d) considers the same population ratio as estimated 
in the experiment [l^, and reproduces quite satisfactorily the observed loss of coherence. 
Guided by panel c), in e) we change Ni from 340 to 400 to probe the opposite phenomenon. 
It turns out that the presence of K enhances V^^-* only at relatively low lattice depths, while 
at large s the effect is again a loss of coherence, albeit less pronounced. This result agrees 
with experiments, where an increase of coherence was never observed 29 1. 

We now turn our attention to another interesting feature of strongly interacting BB mix- 
tures, i.e. the possible occurrence of low energy metastable states characterized by a large 
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number of interfaces, recently discussed in the ideal case of homogeneous lattices, fi/-,- = 0, 

n n ' /.^ ' 

|13l | . The authors of Ref . [13[ observe that the QMC simulations employed to determine the 
ground state of the total Hamiltonian H fail to equilibrate as soon as f/12 > and ascribe 
this behavior to the presence of many low-energy metastable states (where metastable refers 
to the robustness of these configurations against the QMC minimization algorithm, which 
is equipped with nonlocal moves). Being characterized by a large number of interfaces sep- 
arating single-species droplets, these metastable states are dubbed quantum emulsions. The 
relevant energies are found to be linearly dependent on the number of interspecies interfaces. 
One interesting feature of these quantum emulsions is their spontaneous randomness, i.e. the 
fact that the droplets exhibit a disordered spatial arrangement despite the absence of any 
randomness in the Hamiltonian parameters. 

Adopting a self-consistent dynamical search algorithm for the ground state of the homo- 
geneous system in its Gutzwiller form (33|, we find that the BB mixture gets trapped into 
a quantum emulsion state whose energy depends on the number of interfaces, in complete 
analogy with the results obtained in Ref. ISj. This is evident from Fig. [2] illustrating the sit- 



uation on a homogeneous lattice for different values of the hopping to interaction ratios [3^ . 
However, the homogeneous lattice of Ref. is a strongly idealized situation, in which the 
only requirement for the occurrence of quantum emulsions is that Ui 2 be sufficiently larger 
than Ui and t/s Q,Q. 

Moving to the inhomogeneous case typical of actual experimental situations, we find that 
AU = \Ui — U2 \ becomes a further critical parameter for the existence of quantum emulsions. 
This is clearly illustrated in Fig. [3l Panel a) shows the average number of interfaces as a 
function of AU, while the inset is the analogous of the leftmost panel in Fig. [2l Panels 
b) and c) show typical configurations at small and large values of AU, respectively. Note 
that the former is characterized by a significant number of randomly arranged single-species 
droplets. In this case Ji = J2 = O.lf/i, but we obtain similar results also for Ji 7^ J2, 
provided that AU ~ 0. Clearly the number of quantum emulsion states at a given energy 
will be smaller compared to the homogeneous case, owing to the reduced degree of symmetry 
of the system. Indeed, in this case, the energy of each droplet does depend on its position 
in the lattice, due to the local potential contribution. Unlike the one dimensional case, in 
higher dimensions the interface energy of a droplet depends on its size. This fact, along with 
the larger lattice connectivity, is expected to hinder the occurrence of quantum emulsions. 
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FIG. 3: (color online) Quantum emulsions in the presence of a harmonic confinement; a) average 
number of interfaces (over 100 realizations) vs intra-species interaction difference; b) same as Fig. 
[2ja; c) and d) local densities of bosons in two limiting situations. The color key is the same as 
Fig. [Jf) and g). In all cases A''i = = 180, ki = k2 = 10~^f7i, other parameters as in Fig. [2l 



In summary, we have investigated the properties of a strongly interacting bosonic mixture 
loaded into an optical lattice, going beyond the idealized situation of a homogeneous system. 
We considered the inhomogeneities arising from the presence of the harmonic trapping po- 
tential typical of standard experimental setups as well as from the differential gravitational 
sag originated by the difference in the masses of the two bosonic species. We reproduced the 
apparently surprising results of the first experiment involving a BB mixture 12|, providing 
a simple explanation for the observed loss of coherence of ®^Rb in the presence of a small 
fraction of ^^K. Furthermore, our results predict that the opposite phenomenon, i.e. the 
increase of coherence predicted by the "naive argument" proposed in Ref. [l^, is limited 
to sufficiently shallow lattice depths. 

We then investigated the effect of inhomogeneity on the quantum emulsion states formerly 
predicted on homogeneous lattices [13] • In particular, we found that, at variance with the 
homogeneous case, a large value of U12 is not sufficient for their occurrence, and a further 
critical condition is a small value of U2 — Ui. This suggests that the use of Feshbach reso- 
nances could be a crucial ingredient for the experimental observation of quantum emulsions 
in heteronuclear mixtures. An intriguing alternative possibility for the realization of lattice 
BB mixtures with directly built-in conditions Ji = J2 and Ui = U2 could be provided by 
a generalization of the models considered in Ref. |22| , by considering dipolar bosons placed 
on two neighboring lattices with angular relations such that the two sets of atoms interact 
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via a strong interspecies repulsion Ui2- 
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